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Abstract 

We consider a fractional generalization of Hamiltonian and gradient systems. We 
use differential forms and exterior derivatives of fractional orders. We derive fractional 
generalization of Helmholtz conditions for phase space. Examples of fractional gradient 
and Hamiltonian systems are considered. The stationary states for these systems are 
derived. 
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1 Introduction 

Derivatives and integrals of fractional order pi |2] have found many applications in recent 
studies in physics. The interest in fractional analysis has been growing continually during 
the past few years. Fractional analysis has numerous applications: kinetic theories |S1 IH E! , 
statistical mechanics [10111111121, dynamics in complex media jlSl lEl 113 CHI CZI , and many 
others 00111111. 

The theory of derivatives of non-integer order goes back to Leibniz, Liouville, Grun- 
wald, Letnikov and Riemann. In the past few decades many authors have pointed out that 
fractional-order models are more appropriate than integer-order models for various real ma- 
terials. Fractional derivatives provide an excellent instrument for the description of memory 
and hereditary properties of various materials and processes. This is the main advantage of 
fractional derivatives in comparison with classical integer-order models in which such effects 
are, in fact, neglected. The advantages of fractional derivatives become apparent in mod- 
elling mechanical and electrical properties of real materials, as well as in the description of 
rheological properties of rocks, and in many other fields. 

In this paper, we use a fractional generalization of exterior calculus that was suggested in 
|18| ll9j. Fractional generalizations of differential forms and exterior derivatives were defined 
in |TH]- It allows us to consider the fractional generalization of Hamiltonian and gradient 
dynamical systems fK)\ 1^ . The suggested class of fractional gradient and Hamiltonian 
systems is wider than the usual class of gradient and Hamiltonian dynamical systems. The 
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gradient and Hamiltonian systems can be considered as a special case of fractional gradient 
and Hamiltonian systems. 

In section 2, a brief review of gradient systems and exterior calculus is considered to fix 
notation and provide a convenient reference. In section 3, a brief review of fractional (ex- 
terior) calculus will be given to fix notations and provide a convenient reference. In section 
4, a definition of fractional generalization of gradient systems is suggested. In section 5, we 
consider a fractional gradient system that cannot be considered as a gradient system. In 
section 6, we prove that a dynamical system that is defined by the well-known Lorenz equa- 
tions 1211 '^^^ t)e considered as a fractional gradient system. In section 7, a brief review 
of Hamiltonian systems is considered to fix notations and provide a convenient reference. In 
section 8, we consider the fractional generalization of Hamiltonian systems and Helmholtz 
conditions. In section 9, the simple example of fractional Hamiltonian systems is discussed. 
Finally, a short conclusion is given in section 10. 

2 Gradient systems 

In this section, a brief review of gradient systems and exterior calculus PI] is considered to 
fix notations and provide a convenient reference. 

Gradient systems arise in dynamical systems theory [201 1221 1^ • They are described by 
the equation dyi/dt = —gradV{x), where x G -R". In Cartesian coordinates, the gradient 
is given by gradV = SidV/dxi, where x = ejXj. Here and later, we mean the sum on the 
repeated indices i and j from 1 to n. 

Definition 1. A dynamical system that is described by the equations 



is an exact form uj = —dV, where V = V{x) is a continuously differentiable function (0- 
form). 

Here d is the exterior derivative [21]. Let V = V{x) be a real, continuously differentiable 
function on i?". The exterior derivative of the function V is the 1-form dV = dxidV / dxi 
written in a coordinate chart (xi, ...,a;„). 

In mathematics |21j , the concepts of closed form and exact form are defined for differential 
forms by the equation duo = Q for a given form to be a closed form and uo = dh for an exact 
form. It is known, that to be exact is a sufficient condition to be closed. In abstract terms, 
the question of whether this is also a necessary condition is a way of detecting topological 
information, by differential conditions. 





is called a gradient system in if the differential 1-form 
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Let us consider the 1-form (jSj). The formula for the exterior derivative d of differential 
form (j21) is 

where A is the wedge product. Therefore, the condition for u to be closed is 

^ _ ^ = 0. 

In this case, if V{x) is a potential function then dV = dxidV / dxi. The implication from 
'exact' to 'closed' is then a consequence of the symmetry of the second derivatives: 

dxidxj dxjdxi 

If the function V = V{x) is smooth function, then the second derivative commute, and 
equation Q holds. 

The fundamental topological result here is the Poincare lemma. It states that for a con- 
tractible open subset X of -R", any smooth p-form /? defined on X that is closed, is also 
exact, for any integer p > (this has content only when p is at most n). This is not true 
for an open annulus in the plane, for some 1-forms uj that fail to extend smoothly to the 
whole disk, so that some topological condition is necessary. A space X is contractible if 
the identity map on X is homotopic to a constant map. Every contractible space is simply 
connected. A space is simply connected if it is path connected and every loop is homotopic 
to a constant map. 

Proposition 1. If a smooth vector field F = eiFi{x) of system (Op satisfies the relations 

dFi dFj 



dxj dxi 



(4) 



on a contractible open subset X of K^, then the dynamical system (Op is the gradient system 
such that 

dxj _ dV{x) 
dt dxi 

This proposition is a corollary of the Poincare lemma. The Poincare lemma states that 
for a contractible open subset X of -R", any smooth 1-form Q defined on X that is closed, 
is also exact. 

The equations of motion for the gradient system on a contractible open subset X of i?" 
can be represented in the form (jSj). Therefore, the gradient systems can be defined by the 
potential function V = V{x). 

If the exact differential 1-form uj is equal to zero [dV = 0), then we get the equation 

V{x) - C = 0, (6) 

which defines the stationary states of the gradient dynamical system (0). Here C is a 
constant. 
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3 Fractional differential forms 

If the partial derivatives in the definition of the exterior derivative d = dxid/dxi are allowed 
to assume fractional order, a fractional exterior derivative can be defined ^18| by the equation 

= {dx.rT>:^. (7) 

Here, we use the fractional derivative in the Riemann-Liouville form pP that is defined 
by the equation 

where m is the first whole number greater than or equal to a. The initial point of the 
fractional derivative pi is set to zero. The derivative of powers A; of a; is 

where A; > 1, and a > 0. The derivative of a constant C need not be zero 

D^C = ^prf^C. (10) 

For example, the fractional exterior derivative of order a of x^, with the initial point 
taken to be zero and n = 2, is given by 

= (dxi)"D^^a;t + (dx2)"D°^a;^ (11) 

Using equation Q, we get the following relation for the fractional exterior derivative of xf: 

d^x\ = (rfxi)"^^^-^^ + {dx,r 



T{k + l-a) ' ' r(l-a) 

4 Fractional gradient systems 

A fractional generalization of exterior calculus was suggested in [THlll9j . A fractional exterior 
derivative and the fractional differential forms were defined |TB|. It allows us to consider the 
fractional generalization of gradient systems. 

Let us consider a dynamical system that is defined by the equation d^/dt = F, on a 
subset X of R^. In Cartesian coordinates, we can use the following equation 

^ = ^^(^)' (12) 

where i = 1, ..,n, x = e^Xi, and F = eiFi{x). The fractional analog of definition 1 has the 
form. 
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Definition 2. A dynamical system H^) is called a fractional gradient system if the 
fractional differential 1-form 

Uo, = Fi{x){dxiY (13) 

is an exact fractional form uo^ = —d°'V, where V = V{x) is a continuously differentiahle 
function. 

Using the definition of tlie fractional exterior derivative, equation (jl3p can be represented 

as 

uj^ = -d'^V = -{dxi)"'D'^.V. 

Therefore, we have Fi{x) = — D".V^. 

Note that equation (jl3|) is a fractional generalization of equation If a = 1, then 
equation (fT^ leads us to equation (j2)). Obviously, a fractional 1-form uja can be closed when 
the 1-form u = Ui is not closed. The fractional generalization of the Poincare lemma is 
considered in ^^1- Therefore, we have the following proposition. 

Proposition 2. // a smooth vector field F = eiFi{x) on a contractible open subset X of 
R"^ satisfies the relations 

D:^F,-D:^F,=0, (14) 
then the dynamical system / lijjj) is a fractional gradient system such that 

^ = -D^r(x), (15) 
where V{x) is a continuous differentiahle function and D^.l^ = — Fj. 

Proof. This proposition is a corollary of the fractional generalization of Poincare lemma 
[TU] . The Poincare lemma is shown jIHl ^] to be true for the exterior fractional derivative. 

Relations ()14j) are the fractional generalization of relations (j3)). Note that the fractional 
derivative of a constant need not be zero ()10j). Therefore, we see that constants C in the 
equation V{x) = C cannot define a stationary state of the gradient system (jisp . It is easy 
to see that 

T>lV{x) = T)IC = C ^ 0. 

In order to define stationary states of fractional gradient systems, we consider the solutions 
of the system of the equations Y)^y{x) = 0. 

Proposition 3. The stationary states of gradient system Hl^) are defined by the equation 

n m— 1 m— 1 n 

v{x) - lU^r-"^ E - E c......„(n(^o'0 = 0- (16) 

i=l ki=0 k„=0 i=l 
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The Cki...kn o,re constants and m is the first whole number greater than or equal to a. 

Proof. In order to define the stationary states of a fractional gradient system, we consider 
the solution of the equation 

B^Vix) = 0. (17) 

This equation can be solved by using equation (jH}. Let m be the first whole number greater 
than or equal to a; then we have the solution [T1I2] of equation (fTTj) in the form 

m—1 

V{x) = \xi\°''^ak{xi, ...,Xi-i,Xi+i, ...,Xn){xi)'', (18) 

where are functions of the other coordinates. Using equation ()18|) for i = 1, ...,n, we get 
the solution of the system of equation ()17|) in the form p6|) . 

If we consider n = 2 such that x = Xi and y = X2, we have the equations of motion for 
fractional gradient system 

^ = -B»V{x,y), f^=-B-V{x,y). (19) 
The stationary states of this system are defined by the equation 

m— 1 m—1 

r(x,i/) - Ix^r-i J] J] Cfc^xV = 0. 

fc=0 1=0 

The Cki are constants and m is the first whole number greater than or equal to a. 

5 Examples of fractional gradient system 

In this section, we consider a fractional gradient systems that cannot be considered as a 
gradient system. We prove that the class of fractional gradient systems is wider than the 
usual class of gradient dynamical systems. The gradient systems can be considered as a 
special case of fractional gradient systems. 

Example 1. Let us consider the dynamical system that is defined by the equations 

dx ^ dy ^ , . 

-^ = Fx, -T- = Fy, 20 
dt ' dt ^ ' 

where the right hand sides have the form 

= acx^-'' + bx-\ Fy = {ax + b)y-\ (21) 
where a ^ 0. This system cannot be considered as a gradient dynamical system. Using 

oy ox 
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we get that u) = F^dx + Fydy is not closed form 

duj = —ay~^dx A dy. 

Note that the relation (jl4j) in the form 

D°F - D°F =0 

is satisfied for the system (jSTJ, ii a = k and the constant c is defined by 

_ r(l-a) 
r(2-a)' 

Therefore, this system can be considered as a fractional gradient system with the linear 
potential function 

y) = r(l — a){ax + b), 

where a = k. 

Example 2. Let us consider the dynamical system that is defined by equation (I^Uj) with 

= an{n - l)^"-^ + ck{k - l)x''-%\ (22) 

Fy = bm{m - + d{l - l)xV"^ (23) 

where k ^ 1 and 1^1. It is easy to derive that 

dF^ dFy _ k-2„,i-2 



dy dx 



ckl x'-y-'lik - l)y - (/ - l)x] ^ 0, 



and the differential form u = F^dx + Fydy is not closed du ^ 0. Therefore, this system is 
not a gradient dynamical system. Using conditions (jl4p in the form 



y dy'' dx' ' 



we get d°'uo = for a = 2. As the result, we have that this system can be considered as a 
fractional gradient system with the potential function 

V{x, y) = ax" + 6?/" + cx^yK 

In the general case, the fractional gradient system cannot be considered as a gradient 
system. The gradient systems can be considered as a special case of fractional gradient 
systems such that a = 1. 
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6 Lorenz system as a fractional gradient system 

In this section, we prove that dynamical systems that are defined by the well-known Lorenz 
equations [23 fractional gradient system. 



The well-known Lorenz equations j2H[ |21] are defined by 



dx dy dz 



where the right hand sides F^, Fy and F^ have the forms 

Fx = (T{y-x), Fy={r-z)x-y, F^ = xy-hz. 

The parameters cr, r and h can be equal to the following values 

a = 10, h = 8/3, r = 470/19 - 24.74 . 

The dynamical system which is defined by the Lorenz equations cannot be considered as a 
gradient dynamical system. It is easy to see that 

z + a — r, 



dF, 


dFy 


dy 


dx 


OF, 


dF, 


dz 


dx 


dFy 


dF, 


dz 


dy 



-y, 

-2x. 



Therefore, u = F^dx + Fydy + F^dz is not a closed 1-form and we have 

du = —{z + (7 ~ r)dx Ady + ydx Adz + 2xdy A dz. 

For the Lorenz equations, conditions (jl4j) can be satisfied in the form 

DjF.-D^F, = 0, D^F.-D^F, = 0, D^F, - DJF, = 0. 

As the result, we get that the Lorenz system can be considered as a fractional gradient 
dynamical system with potential function 

V{x, y, z) = ^ax^ - ^ayx'^ + ~ ^^^^^ + 

The potential uniquely defines the Lorenz system. Using equation ()16|) . we can get that 
the stationary states of the Lorenz system are defined by the equation 

V{x, y, z) + Coo + C^x + Cyy + C^z + C^yXy + C^^xz + Cy^yz = 0, (25) 
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where Coo, Cx, Cy, Cz Cxy, Cxz, and Cyz are the constants and a = m = 2. The plot of 
stationary states of Lorenz system with the following constants Coo = ^, Cx = Cy = Cz = 
Cxy = Cxz = Cyz = and parameters a = 10, 6 = 3, and r = 25 is shown in figure 1 and 2. 
Note that the Rossler system which is defined by the equations 



dec . . dec dz „ / \ 

— = -{y + z), — =x + 0.2y, — = 0.2+ 
at at at 

can be considered as a fractional gradient system with the potential function 

y{x, y,z) = ^{y + z)x^ - ^xy^ - - ^z^ ^ ~ (^6) 

This potential uniquely defines the Rossler system. The stationary states of the Rossler 
system are defined by equation ^I^. where the potential function is defined by (j^ . The 
plot of stationary states of the Rossler system for the constants Coo = 1, C^,. = C^ = C2 = 
Cxy = Cxz = Cyz = and parameter c = 1 is shown in figures 3 and 4. 

Let us note the interesting qualitative property of surfaces which is difficult to see 
from figures. The surfaces of the stationary states of the Lorenz and Rossler systems separate 
the three dimensional Euclidean space into some number of areas. We have eight areas for 
the Lorenz system, and four areas for the Rossler system. This separation has the interesting 
property for some values of parameters. All regions are connected with each other. Beginning 
movement from one of the areas, it is possible to appear in any other area, not crossing a 
surface. Any two points from different areas can be connected by a curve which does not 
cross a surface. It is difficult to see this property from the figures 1-4. 



7 Hamiltonian systems 

In this section, a brief review of Hamiltonian systems is considered to fix notations and 
provide a convenient reference. 

Let us consider the canonical coordinates (g^, 5",^^, in the phase space R'^"^. 

We consider a dynamical system that is defined by the equations 

The definition of Hamiltonian systems can be realized in the following form fT7\ |2H] ■ 

Definition 3. A dynamical system \21^ on the phase space R^^ , is called a Hamiltonian 
system if the differential 1-form 

(3 = G'dpi - F'dqi, (28) 

is a closed form d(3 = 0, where d is the exterior derivative. A dynamical system is called a 
non- Hamiltonian system if the differential 1-form (3 is nonclosed djS 0. 
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The exterior derivative for the phase space is defined as 

d d 

d = dqi— + dpi—. (29) 
OQi dpi 

Here and later, we mean the sum on the repeated indices i and j from 1 to n. 

Proposition 4. // the right-hand sides of equations satisfy the Helmholtz conditions 



21 . \28^ for the phase space, which have the following forms: 



dpj dpi 



0, (30) 



f^ + I^^O, (31, 

dqi dpj 



dqj dqi 

then the dynamical system \21\} is a Hamiltonian system. 



(32) 



Proof. In the canonical coordinates {q,p), the vector fields that define the system have 
the components {G\F^), which are used in equation ()27|) . Let us consider the 1-form that 
is defined by the equation 

p = G'dpi - F'dqi. 
The exterior derivative for this form can be written by the relation: 

df3 = d{G'dpi) -d{F'dqi). 

It now follows that 

,^ dG' , , dG\ , dF\ , dF' , 
dp = -7:—dqj A dpi + -^r—dpj A dpi - -:T—dq.j A dqi - -:r—dpj A dqi. 
dqj dpj dqj dpj 

Here, A is the wedge product. This equation can be rewritten in an equivalent form as 

,^ fdG^ dF'\ , , 1 fdG^ dG'\ , , 1 fdF' dF^\ , 

Here, we use the skew-symmetry of dqi A dqj and dpi A dpj with respect to the index i and 
j. It is obvious that the conditions (j30|) - (j32|l lead to the equation djS = 0. 

Some of Hamiltonian systems can be defined by the unique function. 

Proposition 5. A dynamical system \21^ on the phase space R^"' is a Hamiltonian 
system that is defined by Hamiltonian H = H{q,p) if the differential 1-form 

(3 = G'dpi - F'dq,, 
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is an exact form (3 = dH, where d is the exterior derivative and H = H{q,p) is a continuous 
differentiable unique function on the phase space. 



Proof. Suppose that the differential 1-form (3, which is defined by equation (j28j) . has 
the form 

^ dH ^ dH ^ 

p = dH = -;r—dpi + -T^dqi. 
opi oqi 

In this case, vector fields (G", F*) can be represented in the form 

dpi dqi 

If if = H{q,p) is a continuous differentiable function, then condition (I^Uj) - are satisfied. 
Using proposition 4, we get that this system is a Hamiltonian system. The equations of 
motion for the Hamiltonian system ()27|) can be written in the form 

dqi _ dH_ dpi^ _ _dl£ ^^^^ 
dt dpi ' dt dqi ' 

which is uniquely defined by the Hamiltonian H. 

If the exact differential 1-form (3 is equal to zero {dH = 0), then the equation 

H{q,p)-C = (34) 
defines the stationary states of the Hamiltonian system (j27p . Here, C is a constant. 



8 Fractional Hamiltonian systems 

Fractional generalization of the differential form ()28|). which is used in definition of the 
Hamiltonian system, can be defined in the following form: 

f3^ = G\dpir-r{dqir. 

Let us consider the canonical coordinates (x^, x", x*^"*"^, x^") = (g^, in 
the phase space i?^" and a dynamical system that is defined by the equations 

The fractional generalization of Hamiltonian systems can be defined by using fractional gen- 
eralization of differential forms ^18j . 

Definition 4. A dynamical system i35]} on the phase space i?^" is called a fractional 
Hamiltonian system if the fractional differential 1-form 

Pa = G^dpiY - F\dq,r, 
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is a dosed fractional form 

d"pa = 0, (36) 

where d°' is the fractional exterior derivative. A dynamical system is called a fractional non- 
Hamiltonian system if the fractional differential 1-form (3^ is a nonclosed fractional form 
d"/3a ^ 0. 

The fractional exterior derivative for the phase space R^"^ is defined as 

d^ = {dqir^l + (dPiT^p.- (37) 

For example, the fractional exterior derivative of order a of q^, with the initial point taken 
to be zero and n = 2, is given by 

d'^q'' = (dg)°D^g'= + (dp)"D;^g^ (38) 

Using equations and ()1U|) . we have the following relation for the fractional exterior deriva- 
tive dnii): 

k Nar(fc + l)g'=-" , ,^ gV" 
T[k + l — a) r(l — a) 

Let us consider a fractional generalization of the Helmholtz conditions. 

Proposition 6. // the right-hand sides of equations / l,?,5]) satisfy the fractional general- 
ization of the Helmholtz conditions in the following form 

^IG^ - D^^G^- = 0, (39) 

D^^G^+D^^F^ = 0, (40) 

D-F^-D-F^=0, (41) 
then dynamical system Vj5\) is a fractional Hamiltonian system. 

Proof. In the canonical coordinates {q,p), the vector fields that define the system have 
the components {G^,F^), which are used in equation (fTTj) . The 1-form j3a is defined by the 
equation 

= G\dpr - F\dq,r. (42) 
The exterior derivative for this form can now be given by the relation 

= d''{G\dpiY) - d''{F\dq,Y). 

Using the rule 

fe=0 
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and the relation 

^{{dxr) = o {k>i), 

we get that 

°° gk 
fc=0 3 

Here, we use 

r(i-a)r(A;+ 1) • 

Therefore, we have 

d"Po. = BlG^dqj f A {dp,r + B'^^G\dp,r A {dp,r- 

-D° A (c/gi)" - B'^^F'idpjy A (dqi)". 

This equation can be rewritten in an equivalent form 

d'^(3. = (t>IG^ + D^^F^) {dqr A {dp,r + \ [^IG^ - D° G^) {dp.T A 

+ \{Dir - BlF=){dq,r ^{dq,r. 

Here, we use the skew symmetry of A. It is obvious that conditions - lead to the 
equation d°'(3oi = 0,i.e., Pa is a closed fractional form. 

Let us define the Hamiltonian for the fractional Hamiltonian systems. 

Proposition 7. A dynamical system on the phase space i?^" is a fractional Hamil- 
tonian system that is defined by the Hamiltonian H = H{q,p) if the fractional differential 
1-form 

l3^ = G\dpir-F\dqi)\ 

is an exact fractional form 

Pa = d^H, (43) 

where is the fractional exterior derivative and H = H{q,p) is a continuous differentiahle 
function on the phase space. 



Proof. Suppose that the fractional differential 1-form Pa-, which is defined by equation 
has the form 

P^ = d^H = (dp.rB^M + {dq^T-DlH. 

In this case, vector fields (G*, F*) can be represented in the form 

G*(g,p) = D^^if, F\q,p) = -T>IH. 
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Therefore, the equations of motion for fractional Hamiltonian systems can be written in the 
form 

^ = B:H, ^ = -BlH. (44) 
dt ' dt ^ ^ 

The fractional differential 1-form (3^ for the fractional Hamiltonian system with Hamil- 
tonian H can be written in the form jSa = d'^H. If the exact fractional differential 1-form (3a 
is equal to zero (d^iJ = 0), then we can get the equation that defines the stationary states 
of the Hamiltonian system. 

Proposition 8. The stationary states of the fractional Hamiltonian system dTTP Q'^e 
defined by the equation 

n m—1 m—1 n 

H{q,p)-\l[q.pr-'^ - E C,,...,j,..,^l[{q,)''{p,y^=0, (45) 

i=l ki=0,li=0 fc„=0,/„=0 1=1 

where Ck^,„k„,ii,...,in '^'^^ constants and m is the first whole number greater than or equal to a. 
Proof. This proposition is a corollary of proposition 3. 



9 Example of fractional Hamiltonian system 

Let us consider a dynamical system in phase space {n = 1) that is defined by the equation 

where the fractional order < a < 1 and the Hamiltonian H[q,p) has the form 

H{q,p) = ap^ + bq^. (47) 

If a = 1, then equation (j46|) describes the linear harmonic oscillator. 
If the exact fractional differential 1-form 

(3a = d'^H = {dpYT>'^H + {dq)'^'D'^H 

is equal to zero ((i"if = 0), then the equation 

H{q,p)-C\qpr~' = 

defines the stationary states of the system Here, C is a constant. If a = 1, we get the 
usual stationary-state equation (j^H). 

Using equation (jTFj) . we get the following equation for stationary states: 

Iqpl^-^iap^ + bq^) = C. (48) 

If a = 1, then we get the equation ap"^ + bq^ = C, which describes the ellipse. 
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10 Conclusion 



Fractional derivatives and integrals ^1 12] have found many applications in recent studies in 
physics. The interest in fractional analysis has been growing continually during the past few 
years [3-17]. Using the fractional derivatives and fractional differential forms, we consider 
the fractional generalization of gradient and Hamiltonian systems. In the general case, the 
fractional gradient and Hamiltonian systems cannot be considered as gradient and Hamil- 
tonian systems. The class of fractional gradient and Hamiltonian systems is wider than the 
usual class of gradient and Hamiltonian dynamical systems. The gradient and Hamiltonian 
systems can be considered as aspecial case of fractional gradient and Hamiltonian systems. 
Therefore, it is possible to generalize the application of catastrophe and bifurcation theory 
from gradient to a wider class of fractional gradient dynamical systems. Note that quan- 
tization of the fractional Hamiltonian systems can be realized by the method suggested in 

[2niEniEIlE21. 
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